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ABSTRACT 



Context. Stellar mergers are expected to take place in numerous circumstences in the evolution of stellar systems. In particular, they 
are considered as a plausible origin of stellar eruptions of the V838 Mon type. V1309 Sco is the most recent eruption of this type in 
our Galaxy. The object was discovered in September 2008. 
Aims. Our aim is to investigate the nature of VI 309 Sco. 

Methods. V1309 Sco has been photometrically observed in course of the OGLE project since August 2001. We analyse these obser- 
vations in different ways. In particular, periodogram analyses were done to investigate the nature of the observed short-term variability 
of the progenitor. 

Results. We find that the progenitor of VI 309 Sco was a contact binary with an orbital period of ~ 1.4 day. This period was decreasing 
with time. The light curve of the binary was also evolving, indicating that the system evolved towards its merger. The violent phase of 
the merger, marked by the systematic brightenning of the object, began in March 2008, i.e. half a year before the outburst discovery. 
We also investigate the observations of VI 309 Sco during the outburst and the decline and show that they can be fully accounted for 
within the merger hypothesis. 

Conclusions. For the first time in the literature we show from direct observations that contact binaries indeed end up by merging into 
a single object, as was suggested in numerous theoretical studies of these systems. Our study also shows that stellar mergers indeed 
result in eruptions of the V838 Mon type. 

Key words, stars: individual: V1309 Sco - stars: binaries: general - stars: peculiar 



1. Introduction 

Stellar mergers have for a long time been recognized to play 
an important role in the evolution of stellar systems. High stel- 
lar densities in globular clusters can lead often to collisions and 
mergers of stars (Leonard 1989). In this way the origin of blue 
strugglers can be explained. In dense cores of young clusters 
multiple mergers of protostars have been suggested as a way of 
formation of the most massive stars (iBonnell et al.lll998l) . Some 
binary stars, in particular contact binaries, are suggested to end 
their evolution as stellar mergers (Robertson & Egaleton 1977). 

The powerful outburst of V838 Mon in 2002 ^Munari et al. 
2OO2I) . accompanied by a spectacular light echo (Bond et aL 



which erupted in 2008 (iMason et al. l l2oToh . As extragalactic 
erupti ons of this kind o ne can mention M31 RV (eruption i n 
1989, 'Mould et al."l990'), M85 OT2006 (Kulk arni et aT]l2007h . 
and NGC300 OT2008 (Bergeretal. 2009). Several interpreta- 
tions of the er uptions were proposed . They include an u nusual 
classical nova (llben & T utuW' 19921 IShara et al.ll2010h . a late 



20031) . raised interest in a class of stellar eruptions named "red 
novae", "optical transients" or "V838 Mon type eruptions". 
These objects, which typically reach a maximum luminosity 
of ~ 10* Lq, evolve to low effective temperatures and de- 
cline as very cool (super)giants. Apart form V838 Mon, in 
our Galaxy the class also includes V4332 S gr, whose outburst 
was observed in 1994 (iMartini et al.l Il999h . and V1309 Sco, 
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* Based on observations obtained with the 1.3-m Warsaw tele- 
scope at the Las Campanas Observatory of the Carnegie Institution of 
Washington. 



He-shell flash (lLawloij l2005), or a thermonuclear shell flash in 
an evolved massive star (Munari et al.ll2005l) . iTvIenda & Sokeil 
(2006) presented numerous arguments against these mecha- 
nisms. They showed that all the main observational characteris- 
tics of the V838 Mon-type eruptions can be consistently under- 
stood as resulting from stellar collis i ons a nd mergers, as orig- 
inally proposed i n Soker & Tvlendal (12003). For these reasons 
ISoker & Tvlendal (12007) proposed to call these t ype of eruptions 
mergebursts. Recently, iKashi et aP (|2010) and iKashi & Soked 
have suggested that some of the V838 Mon-type erup- 
tions are of the same nature as the eruptions of luminous blue 
variables, and that they can be powered by mass-transfer events 
in binary systems. 

In the present paper, we show that the recent red nova, 
VI 309 Sco, is a Rosetta stone in the studies of the nature of 
the V838 Mon type eruptions. Archive photometric data col- 
lected for the object in the OGLE project during about six years 
before the outburst allow us to conclude that the progenitor of 
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Fig. 1. Light curve of VI 309 Sco from the OGLE-III and OGLE-IV projects: / magnitude versus time of observations in Julian 
Dates. Time in years is marked on top of the figure. At maximum the object attained / - 6.8. 



VI 309 Sco was a contact binary. The system quickly evolved 
towards its merger, which resulted in the eruption observed in 
2008. 

2. V1309 Scorpii 

VI 309 Sco, also known as Nova Sco 2 008, wa s disco vered 
on 2.5 September 2008 (JD 2454712) (iNakanol l2008h . The 
subsequent evolution, however, showed that as pointed out in 
[Mason et al. (2010), this was not a typical classical nova. 

Ea rly spectroscopy revealed an F-type giant (Ru dv et alj 
I2008ah .' On a time scale of a month the object evolved to K- 



and early M-types (iMason et al bolOtlRudv et al.ll2008bh . Eight 
months after t he discovery it was observed as a late M-type giant 
([Mason et al ] [2oioh . 

As described in iMason et al.l (1201 O l). the object developed a 
complex and rapidly evolving line spectrum of neutral and singly 
ionized elements. In early epochs, absorption features dominated 
the spectrum, emission components developed later and their 
intensities quickly increased with time. The strongest emission 
lines were those of the hydrogen Balmer series. 

The emission components had FWHM of ~150 kms ' and 
broader wings , which in the case of Ha extended even beyond 
1000 kms ' dMason et al.l l2010l) . Narrow absorption compo- 
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nents were superimposed on the emission components, so that 
the line profiles mimicked those of P-Cyg type or inverse P-Cyg 
in some cases. Mason et al. (2010) interpret these line profiles 
as produced in an expanding shell that is denser in the equatorial 
plane. 

In late epochs, when the object evolved to M-type, absorp- 
tion features of TiO, VO, CO, and H2O appeared in the spectrum 
(Mason et al. 2010; Rudy et al. 2008b). 

VI 309 Sco thus shares the principal characteristic of the 
V838 Mon type eruptions, i.e evolution to very low effective 
tempratures after maxim um brightness and during the decline 
(iTvlenda & Sokeiil2006h . Other common features of VI 309 Sco 
and the V838 Mon type eruptions include: outburst time scale 
of the order of months, outburst amplitude (maximum minus 
progenitor) of 7 - 10 magnitudes, complete lack of any high- 
ionization features (coronal lines, in particular), expansion ve- 
locities of a few hundred km s"' (instead of a few thousands as 
in classical novae), and oxide bands observed in later epochs, 
which imply that oxygen-rich (C/O < 1 ) matter was involved in 
the eruption. 

3. Observations 

Owing to the position of VI 309 Sco close to the Galactic centre 
(1 = 359?8, b = -3^1), the object appears to be situated within a 
field m onitored in the OGLE-III and OGLE-IV projects ( Udalski 
I2003IP1 . As a result VI 309 Sco was observed on numerous occa- 
sions since August 2001. Altogether more than 2000 measure- 
ments were obtained predominantly in the / Cousins photomet- 
ric band. Among them, ~1340 observations were made before 
the discovery of the object as a nova in September 2008. A few 
observations were also made in the V band (seven before the 
discovery). The data were reduced and calibrated using standard 
OGLE procedures (Udalski et al. 2008). A typical precision of 
the measurements was 0.01 magnitude. 

4. Results 

The entire measurements of V 1 309 Sco derived from the OGLE- 
III and IV surveys in the / photometric band are displayed in 
Fig.[T] The gaps in the data are owing to conjunctions of the ob- 
ject with the Sun. Apart from 2001 and 2009 most of the data 
were obtained between February and October of each year. In 
2001 the OGLE-III project was just starting to operate, hence 
only 11 measurements were obtained in this year (they were 
omitted from the analysis discribed below). In the period 2002- 
2008, from 52 (2002, 2003) to 367 (2006) observations were 
made each year Near the maximum of the 2008 eruption, i.e. 
when the object was brighter than / ^ 11, its image was over- 
exposed in the OGLE frames, hence there are no data for this 
period. Near maximum brightness, the object attained / ^ 6.8, 
according to the data gathered in AAVSCB In May 2009 OGLE- 
III phase ended and threfore only 64 data points were collected 
during the 2009 observing season. The OGLE project resumed 
regular observations of the Galactic centre in March 2010 with a 
new instrumental setup, namely a 32 chip mosaic camera, which 
started the OGLE-IV phase. In 2010 655 measurements were 
obtained. 

As can be seen from Fig. [1] the progenitor of VI 309 Sco 
was initially slowly increasing in brightness on a time scale of 
years and reached a local maximum in April 2007. Subsequently, 

' http://ogle.astrouw.edu.pl 

^ available at http://www.aavso.org/| 



the star faded by ~1 magnitude during a year In March 2008, 
the main eruption started, which led to the object discovery, as 
Nova Sco 2008, half a year later. Remarkable is the smooth, 
roughly exponential rise in brightness. During this event the ob- 
ject brightened by ~10 magnitudes, i.e. by a factor of 10'*. An 
analysis of the data for the progenitor (2002-2007) is presented 
in Sect. |5] while the outburst and the decline of the object are 
discussed in Sect.|6l 

5. The progenitor 

5.1. The data 

A short-term variability, resulting in a ~0.5 magnitude scatter 
of the points in Fig. [1] is the most remarkable and interesting 
feature of the VI 309 Sco progenitor This variability is strictly 
periodic. To show this, we used a method that employs periodic 
orthogonal polynomials to fit the observations and an analysis 
of variance to evaluate the quality of the fit, as described in 
ISchwarzenberg-Czernvl d 1 99 6). This method is particularly suit- 
able for analysing unevenly sampled data of non-sinusoidal peri- 
odic variations, such as eclipses in binary systems. The resulting 
periodograms obtained for the particular observational seasons 
are presented in Appendix lAl They show that the principal pe- 
riodicity in the variations observed in seasons 2002 - 2007 cor- 
responds to a period of ~ 1 .4 day. However, the derived period 
was not constant, but slowly decreasing with time, as presented 
in Fig.|2] It decreased by 1.2% during the pre-outburst observa- 
tions. 

When deriving the period plotted in Fig. |2l we devided the 
data sequence available for a given season into subsamples, each 
containing -50 data points, and constructed periodograms for 
each subsample separately. In this way we can see that during 
later seasons, particularly in 2006 and 2007, the period was sig- 
nificantly varying on a time scale of months. 

We find that the time evolution of the observed period can 
be fitted by an exponential formula. A result of a least-squares 
fit (inverse errors used as weights) is shown with the full line in 
Fig.|2] The obtained formula is 

/ 15 29\ 

P = 1.4456 exp ^^— , (1) 

where P is the period in days, f is the time of observations in 
Julian Dates, and fo = 2455233.5. 

Using the above period fit we folded the observations and 
obtained light curves of the object in particular seasons. The re- 
sults are displayed in Fig. [3] In 2002 - 2006 the light curves were 
obtained from all available data for a particular seasons. The re- 
sults are shown in the upper part of Fig. [3] In 2007 the light curve 
significantly evolved during the season and therefore we plotted 
the results for the subsamples (the same as when deriving the pe- 
riod) of 2007 season in the lower part of Fig. |3](time goes from 
subsample a to e in the figure). 

As can be seen from Fig. |3] the light curve in 2002 - 2006 
displays two maxima and two minima. They are practically 
equally spaced in the phase and have similar shapes in the early 
seasons. This is the reason why the periodograms for these sea- 
sons (see Appendix lAli show a strong peak at a frequency twice 
as high (a period twice as short) as the main peak. In other words, 
the observations could have been interpreted with a period of 
~0.7 day and a light curve with one maximum and minimum. 
The periodograms show however already in the 2002 season that 
the period that is twice as long (~1.4 day) better reproduces the 
observations. The reason is obvious in the later seasons, when 
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Fig. 2. Evolution of the period of the photometric variations of 
the VI 309 Sco progenitor. The Hne shows a least-squares fit of 
an exponential formula to the data (see text and Eq.[T]l. 

the first maximum (at phase 0.25 in Fig.O in the light curve be- 
comes increasingly stronger than the second one. As a result the 
peak corresponding to the 0.7 day period deacreases and practi- 
cally disappears in 2007. 

5.2. Interpretation 

In principle one can consider three possible interpretations of the 
observed light variability, i.e. stellar pulsation, single-star rota- 
tion, and an eclipsing binary system. 

5.2.1. Pulsation 

The light curve in the early seasons may imply a stellar pulsation 
with a period of ~0.7 day and an amplitude of ~0. 15 magnitude. 
This interpretation, however, encounters severe problems when 
it is used to explain the observed evolution of VI 309 Sco. As we 
show below, the progenitor was probably a K-type star There 
is no class of pulsating stars that could be reconciled with these 
characteristics, i.e a K-type star pulsating with the above period 
and amplitude. Moreover, to explain the observed evolution of 
the light curve displayed in Fig. [3] one would have to postulate 
a switch of the star pulsation to a period exactly twice as long 
on a time scale of a few years, together with a gradual shorten- 
ing of the period. This would be very difficult, if not impossi- 
ble, to understand within the present theory of stellar pulsation. 
Finally, there is no physical mechanism involving or resulting 
from a pulsation instability, which could explain the powerful 
2008 eruption, when the object brightened by a factor of 10'*. 

5.2.2. Rotation of a single star 

The observed light curve could have been explained by a sin- 
gle K-type star rotating with a period of ~ 1 .4 day and with two 
spots (or rather two groups of spots) on its surface in early sea- 
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Fig. 3. Light curves obtained from folding the data with the pe- 
riod described by Eq. ([TJ. Upper part: seasons 2002 - 2006. 
Lower part: season 2007 devided into five subsamples (time goes 
from a to e). The zero point of the magnitude (ordinate) scale is 
arbitrary. 



sons. The two spots would be replaced by one spot in 2007. The 
object would thus belong to the FK Com class of rapidly rotating 
giants (e.g. Bopp & Stencel 1981). There are, however, several 
differences between the progenitor of VI 309 Sco and the FK 
Com stars. The 1.4 day period is short, as for the FK Com stars. 
The very good phase stability of the observed light curve, shown 
in Fig. [3] implies that the spot(s) would have had to keep the 
same position(s) on the star surface over a time span of a few 
years. This seems to be very improbable given differential rota- 
tion and meridional circulation, which are expected to be sub- 
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stantial for a rapidly rotating star. Indeed, for the FK Com stars, 
the spots are observed to migrate and change their position on 
much shorter time scales (Korhonen et al. 2007). The systematic 
decrease of the rotational period is also difficult to explain for a 
single star, and is not observed in the FK Com stars. Finally, no 
eruption such as that of VI 309 Sco in 2008 was observed for the 
FK Com stars. Indeed, there is no known mechanism that could 
produce such a huge eruption in a single giant, even if it were 
fast rotating. 

5.2.3. A contact binary evolving to its merger 

We are thus left with an eclipsing binary system. As we show 
below, this possibility allows us to explain all principal charac- 
teristics of the observed evolution of the VI 309 Sco progenitor, 
as well as the 2008 outburst. 

The shape of the light curve, especially in early seasons 
(two rounded maxima of comparable brightness and two equally 
spaced minima), implies that the progenitor of VI 309 Sco was 
a contact binary. The orbital period of ~1.4 day does not al- 
low us to classify the object as a W UMa-type binary, because 
the clasical W UMa stars have periods < 1 day. Nevertheless, 
contact binaries with periods beyond 1 day are also observed 
dPaczvriski et al.l [2006). The exponentially decreasing period 
(Fig. |2]i can be interpreted as resulting from an unstable phase 
of evolution of the system, which leads to the shrinkage of the 
binary orbit and finally to the merger of the system. Such a 
situation was predicted in theoretical studies (IWebbinklll976l: 
[Robertson & Eggleton 1977; Rasio 1995) and it is now com- 
monly believed that the W UMa binaries end their evolution by 
merging into a single star Dissipation of the orbital energy in 
the initial, violent phase of the merger dSoker & Tvlendall2006h 
resulted in the VI 309 Sco eruption observed in 2008. 

One of the possible ways that can lead a binary system to 
merge is the so-called Darwin instability. This happens when the 
spin angular momentum of the system is more than a third of the 
orbital angular momentum. As a result, tidal interactions in the 
system cannot maintain the primary component in synchroniza- 
tion anymore. The orbital angular velocity is higher than the pri- 
mary's angular velocity, the tidal forces increase their action and 
rapidly transport angular momentum from the orbital motion to 
the primary's rotation. For contact binaries this takes place when 
the binary mass ratio q = M2/M1 ^ 0.1 (Rasio" 19951) . 

Another possibility appe ars to be if a binary system enters a 
deep contact, as discussed in lWebbinkl(ll976h . This can happen, 
for instance, if the primary attempts to cross the Hertzsprung 
gap. The system then starts loosing mass and angular momen- 
tum through the outer Lagrangian point, L2- This shrinks the bi- 
nary, which further deepens the contact and increases mass and 
angular momentum loss. In addition the system starts orbiting 
faster than the components rotate. As a result, similar as for the 
Darwin instability, the tidal forces transport angular momentum 
from the orbit to the components' spins, which further acceler- 
ates the orbit shrinkage. 

A maximum in the light curve of a typical contact binary is 
observed when we look at the system more or less perpendicu- 
larly to the line that joins the two components. As a result, we 
observe two maxima during each orbital revolution. The max- 
ima are of similar brightness, because the system looks similar 
from both sides. The situation changes when one of the above 
instabilities sets in. The secondary starts orbiting faster than the 
primary's envelope rotates. The stars are in contact, which means 
that the difference in velocities is partly dissipated near the con- 
tact between the components. This should lead to the formation 
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Fig. 4. V (open points) and / (asterisks) measurments obtained in 
JD 2453880 - 2453910 (season 2006) and folded with the period 
described by Eq. ([TJ. The zero point of the magnitude (ordinate) 
scale is arbitrary. 



of a brighter (hotter) region that is visible when we look at the 
leading side of the secondary. The system is then brighter than 
when we look from the opposite direction. The maxima in the 
light curve start to differ. This is what we observe in the case 
of the VI 309 Sco progenitor As can be seen from Fig. |3] the 
first maximum (at phase ~0.25) becomes progressively stronger 
than the second one (phase -0.75) from 2002 to 2006. In 2007, 
the second maximum disappears, and at the end of this season 
we observe a light curve with only one maximum and one min- 
imum. Apparently the system evolved to a more spherical con- 
figuration with a bright (hot) spot covering a large fraction of the 
system's surface along the orbital plane. 

5.3. Basic parameters of the progenitor 

5.3.1. Interstellar reddening, spectral type, and effective 
temperature 

The effective temperature of the progenitor can be derived from 
the V - I colour. In 2006 a few measurements were obtained in 
the V band. They are displayed in Fig.|4]together with the / data 
obtained in the same time period (JD 2453880 - 2453910). 

In order to determine the effective temeprature from a 
colo ur, an estimat e of th e interstellar reddening is also neces- 
sary. iMason et"al ] (120101) obtained Eb-v - 0.55 from inter- 
stellar lines observed in their spectra of VI 309 Sco in out- 
burst. We constrained spectral ty pes of the object in outburst 
from the UVES spectra of iMason et al.l (1201 Ol) obtained in their 
epochs 1-3 (downloaded from the VLT archive). Comparing 
them with the BVRI photometry of the object obtained by the 
AAVSO team near the above spectral observations, we obtained 
0.7 $ Eb-v ~ 0.9. The OGLE data also allowed us do derive 
a V - I colour, equal to 5.27, 223 days after the outburst dis- 
covery date. This happened to be in the epoch of the SOAR 



6 



Tylendaetal: V1309 Sco 



spectroscopy of iMason et al ] (1201 Oh . from which these authors 
derived a spectral type of M6-7. To reconcile the colour with 
the spectral type, Eb-v > 0.7 is required. An upper limit to the 
extinction can be obtained from the Galactic Dust Extinction 
Servic^E which estimates Eb-v from the 100 yum dust emis- 
sion mapped by IRAS and COBE/DIRBE. For the position of 
VI 309 Sco, we obtain Eb-v ~ 1.26 ± 0.11. Below we assume 
that V1309 Sco is reddened with Eb-v - 0.8. 

As can be seen from Fig.|4l two V measurements were taken 
close to the first maximum in the light curve (phase ~0.26), two 
others near the following minimum (phase ~0.55). The resulting 
colour, V -I, was 2.09 and 2.21 for the two phases, respectively. 
The mean value of the observed colour, i.e. V - I - 2.15, cor- 
rected for the reddening and compared with the standard colours 
of giants (luminosity class III) results in a spectral class Kl-2 
and an effective temperature Teff - 4500 K for the progenitor of 
VI 309 Sco. 

The difference in the colour, as mentiond above, implies that 
at the maximum the object was on average ~200 K hotter than 
in the minimum. This agrees with the interpretation of the evo- 
lution of the hght curve given in Sect. 15.2.31 

5.3.2. Luminosity and distance 

The hypothesis that the progenitor of VI 309 Sco was a con- 
tact binary allows us to evaluate other parameters of the system. 
Adopting a total mass of the binary, (Mi + M2), of 1.0 - 3.0 Mq 
(observed range of masses of W UMa binaries) and an orbital 
period of 1.43 day, one obtains from the third Kepler law that 
the separation of the components. A, is (3.7 - 5.4) x 10" cm (5.4 
-7.7Ro). 

Taking the mass ratio, q, within the range observed in the 
W UMa systems, i.e. 0.07 - 1 . 0, we can derive the radii of the 
Roche lobes of the components (Egg letonll983b . Assuming that 
both components fill their Roche lobes a total effective surface of 
the system can be obtained. This, along with the above estimate 
of Teff results in a luminosity of the system of (1.15 - 3.3) x 
lO^'* ergs-i (3.0- 8.6 Lo). 

Taking into account bolometric correction, interstellar red- 
dening, and the mean observed brightness of the VI 309 Sco pro- 
genitor in 2006, V = 18.22, we derive a distance to the object of 
3.0 + 0.7 k pc. If the outer Roche lob es (dimensions can be calcu- 
lated from lYakut & Eggletonll2005l) are assumed to be filled by 
the components, the luminosity increases by ~35%, while the 
distance becomes 3.5 + 0.7 kpc. The distance estimate is almost 
independent of the reddening value. 

The above estimates of the luminosity and effective temper- 
ature are consistent with a ~1.0 star at the beginning of the 
red giant branch when they are compared with theoretical tracks 
of stellar evolution (e.g. ,Girardi et al.,.2000,) . 



6. The outburst 

6.1. The data 

Figure |5] (upper part) displays the rise of VI 309 Sco during its 
eruption in 2008 (open points near maximum show the / results 
taken from AAVSCO. This phase started in Mai'ch 2008, and it 
took the object about six months to reach the maximum bright- 
ness. The rise was remarkably smooth for an eruption that led 
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Fig. 5. / light curve of VI 309 Sco during its rise to maximum in 
2008 (upper part). Full points: data from OGLE 111, open points: 
data from AAVSO. The line shows a fit of an exponential for- 
mula (Eq.|2]i. The time scale used in the formula is plotted in the 
lower part of the figure. 

the object to brighten by a factor of ~ 10"*. We found that the 
light curve can be fitted by an exponential formula, i.e. 

/ = -2.5 log(F/), Fj^Fo+ Fi exp(f/T(f)), (2) 

where Fo - l.6x 10"^, Fi = 2.1 x 10"*^ are fluxes in the Vega 
units, f = JD - 2454550, and T{t) is a time scale, allowed to 
evolve with time. The time evolution of t, resulting in the fit 
shown with the full line in the upper part of Fig. |5] is displayed 
in the lower part of the figure. These results show that during 
the initial ~5 months of the rise, i.e. up to JD ^ 2454690, the 
time scale was constant and equal to -27 days. During this phase 
the object brightened by ~3.5 magnitudes. Then the time scale 
quickly decreased and the object brightened by ~6.5 magnitudes 
during a month. 
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Table 1. Basic parametres of VI 309 Sco in outburst and decline 
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Fig. 6. / light curve of VI 309 Sco during the decline in 2010. 



The / light curve of VI 309 Sco during the decline observed 
by OGLE-IV in 2010 is presented in Fig.|6] Evidently the decline 
was relatively smooth at the begining of the presented phase. A 
similar behaviour (smooth evolution) was also observed in 2009 
(see Fig. [TJ. However, as the decline continues in 2010, oscil- 
lations of a few hundredths magnitude on a time scale of hours 
appear A periodogram analysis shows that there is no significant 
periodicity in these light variations. 

In a few cases multiband photometric observations are avail- 
able during the outburst and decline. They allow us to estimate 
the spectral type and effective temperature of the object. This 
is the case for seven dates in September 2008, when multi- 
band photometry was made by the AAVS O team. On 4 and 
16 October 2008 the object was observed bv lRudv et al] (l2008bh 
using the Infrared Telescope Facility, and JHK magnitudes were 
derived. These results can be combined with BV AAVSO mea- 
surements and an / OGLE magnitude obtained in the same time 
period. On 13 April 2009, as well as in a few cases in March- 
April 2010 and October 2010, V photometric measurements (to- 
gether with /) were obtained by OGLE. On 28 August 2010 the 
object was observed by one of us (M.H.) with the SAAO 1.0 
m. Elizabeth Telescope with the Bessell BVRI filters. The ob- 
servations were reduced using standard procedures in the IRAF 
package. The photometry of the star was obtained with the 
DAOPHOT package implemented in IRAF. The resulting mag- 
nitudes are B = 22.21+0.10, y = 20.42+0.13,/; = 18.84+0.10, 
and / = 16.67 ± 0.10. The latter value agrees very well with an 
OGLE result (/ = 16.69 ± 0.01) obtained on the same date. 

6.2. Observed evolution 



Taking our estimates of the distance (3.0 kpc, see Sect. 15.3.2b 
and reddening {Eb-v - 0.8, see Sect. 15.3.1b . we can derive ba- 
sic parameters of VI 309 Sco during the outburst from multi- 
band photometry with the same method as in Tylenda (2005). 
The resulting values of the eff'ective temperature, radius, and lu- 
minosity are presented in Table [1] and Fig. [T] They show that 
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the evolution of VI 309 Sco during outburst and decline wa s re- 
ally of the s ame sort as those of V838 Mon (lTvlendall2005h and 
V4332 Sgr (' Tvlenda et al.l l2005V In all these cases the main de- 
cline in luminosity was accompanied by a decline in the effective 
temperature. In later phases of the decline the objects resume a 
slow increase in the effective temperature. There is no doubt that 
V1309 Sco was not a classical nova. We know from the present 
study that its eruption resulted from a merger of a contact binary. 
Below we show that the energy budget of the eruption can also 
be well accounted for by dissipation of the orbital energy of the 
progenitor. 

6.3. Merger-powered outburst - mergerburst 

The orbital energy and angular momentum of a binary system 
can be calculated from 



^orb 



and 



GM1M2 G (Ml + M2f q 



2 A 



2A 



= ir. 77- MiM2^ (G (Ml + M2? A) 

\Mi+ M2I ^ ' 



V2 q 



(3) 



-.(4) 



Taking the values assumed and obtained in Sect. 15.3.21 for the 
VI 309 Sco progenitor, one obtains iSoi-b - (0-2 - 5.6) 10"*^ ergs 
and 7o,-b = (0.9 - 22.) 10^' g cm^ s"'. 

As derived in Sect. 15. II the orbital period decreased by 1.2% 
during ~6 years. Thus the system contracted by 0.8% and its 
orbital energy decreased by ~ 10"*^ ergs during this time period, 
which gives a mean rate of energy dissipation of ~ 10^*^ erg s"'. 
Already a small portion from this can account for the observed 
brightening of the object between 2002 and the beginning of the 
2007 season (see Fig.[T]l. 

During the same time period the system also lost 0.4% of 
its orbital angular momentum. If the system started shrinking 
because of the Darwin instability, the orbital angular momentum 
loss was primarily absorbed by the spinning-up primary. If, how- 
ever, the system entered a deep contact resulting from evolution- 
ary changes of the components, mass loss through the L2 point 
would be the main source of the orbital angular momentum loss. 
Assuming that the outflowing mass carries out a specific angular 
momentum of ^r^^, where Q is the angular rotational velocity 
of the system and r^,, is a distance of L2 from the mass cen- 
tre, we can estimate that the system lost (0.16 - 2.0)10""' Mq in 
2002-2007. This is very much an upper limit for the mass loss. 
First, judging from the K spectral type, the system might have 
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Fig. 7. Evolution of VI 309 during outburst and decline (see 
Table [TJ. The abscissa displays a logarithm of time in days 
counted from the date of the discovery. 



had magnetically active components. Then the outflowing matter 
might have kept corotation to radii larger than r/,, and thus have 
carried out a specific angular momentum larger than assumed in 
the above estimate. Second, matter already lost through L2 may 
have formed an excretion disc around the system. Tidal interac- 



tions between the binary system and the disc can transfer angular 
momentum directly from the system to the disc. Third, a part of 
the orbital angular momentum probably went to accelerating the 
spins of the components. 

An excretion disc is expected to be formed even if the orbital 
shrinkage was initiated by the Darwin instability. Decreasing 
separation between the components results in a deeper contact, 
and mass loss through L2 must finally occur 

Apparently we see in the light curve of the VI 309 Sco 
progenitor a signature of an excretion disc formation. We ob- 
served the system near the equatorial plane (eclipsing binary). 
Therefore, when the disc became sufficiently massive and thick, 
it could result in dimming the object for the observer. This is a 
likely explanation of the fading of the VI 309 Sco progenitor ob- 
served during a year before the beginning of the 2008 eruption 
(see Fig.[T]). 

The progressive shortening of the orbital period finally led to 
engulfing the secondary by the primary's envelope. This proba- 
bly took place at some point in February 2008, when the signs 
of the binary motion disappeared from the light curve. The sec- 
ondary, or rather its core, now spiralling in the common enve- 
lope, started to release the orbital energy and angular momen- 
tum at an increasing rate. The result was the gradual and rela- 
tively gentle brightening of the object, doubling the brightness 
every ~19 days (Fig.|5]). This phase lasted until ~20 August (JD 
^2454700), when the eruption abruptly accelarated and the ob- 
ject brightened by a factor of ~300 during ~10 days. Perhaps 
this was a signature of a final disruption of the secondary's core 
deeply in the envelope. 

In the merger process, especially during its initial, relatively 
gentle phases until about August 2008, mass loss probably oc- 
cured mainly in directions close to the orbital plane of the pro- 
genitor. As a result, it is likely that an extended disc-like enve- 
lope was then formed, where a significant portion of the angu- 
lar momentum of the progenitor was probably stored. The main 
eruption in August 2008, partly blocked by the e nvelope, was 
then U kely to occur mainly along the orbital axis. iMason et al.l 
ilOlO) interpret the line profiles observed during the outburst, 
especially those of the Balmer series, as produced in a partially 
collimated outflow and a slowly expanding shell that is denser in 
the equatorial plane. 

As can be seen from Table [1] VI 309 Sco attained a max- 
imum luminosity of ~ 3 10* Lq. The OGLE data show that 
the period when VI 309 Sco was brighter than 7^11 mag- 
nitude lasted 40 days. Both the maximum luminosity and the 
time scale of the outburst are close to the theoretical expecta- 
tions of a solar-mass s tar merging with a low-mass companion 
(ISok er & Tvlenda 2006). During ~30 days, when the luminos- 
ity of VI 309 Sco was ^10"* Lq, the object radiated energy of 
~ 3 10'*'* ergs. This is ^1% of the energy available in the bi- 
nary progenitor (see above). Of course a considerable energy 
was also lost during the six-months rise, as well as in the de- 
cline after maximum. Certainly a significant amount was carried 
out in mass loss. Finally, an energy was also stored in the inflated 
remnant. Estimates made for V838 Mon show that the total en- 
ergy involved in the outburst can be a factor of 10 - 20 higher 
than the energy observed in radiation (Tvlenda & Soker 200<^. 
Even this case can be accounted for by the available orbital en- 
ergy. 

Our SAAO photometry performed in August 2010, as well 
as the OGLE V and I measurements obtained in 2010 (see 
Sect. 16. II ). i.e about two years after the outburst maximum, show 
that VI 309 Sco reached a luminosity (see Table [1) comparable 
with its preoutburst value (see Sect. I5.3.2] |. The object was then 
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significantly cooler however than before the outburst (early M- 
type spectrum versus Kl-2). The latter agrees with what was ob- 
served in V838 Mon and V4332 S gr, disp laying M-ty pe sp ectra 
several years after their outbursts dTvlenda 2005; Tvlenda et al.l 
l2005h . However, the drop in luminosity of VI 309 Sco was un- 
expectedly deep. Both V838 Mon and V4332 Sgr at present 
(many years after outburst) remain significantly more luminous 
than their progenitors. Most likely the remnant of VI 309 Sco 
contracting after the outburst partly disappeared for us behind 
the disc-like envelope. Moving small-scale blobs of the enve- 
lope matter, absorbing and scattering the light from the central 
object, could have been responsible for the short-term variabil- 
ity shown in Fig. |6] The situation thus looks to be similar to 
that of the V4332 Sgr remnant, where th e central object is most 
likely hidden in a n opaque dusty disc (iKaminski et alj 1201 Ot 
iKaminski & Tvlenda 2011 ) . 

7. Conclusions 

The principal conclusion of our study is that all observed prop- 
erties of VI 309 Sco, i.e. the light curve of the progenitor during 
six years before the 2008 eruption, as well as the outburst itself, 
can be consistently explained by a merger of a contact binary. 
This is the first case of a direct observational evidence show- 
ing that the contact binary systems indeed end their evolution by 
merging into single objects, as predicted in numerous theoretical 
studies of these systems. 

Our study also provides a conclusive evidence in favour 
of the hypothesis that the V838 Mon-type eruptions (red no- 
vae) result frorn stellar me rgers, as original l y proposed in 
ISoker & Tvleiidal (12003) and iTvlenda & Soke3 0006.) . In pai-- 
ticular, long- and short-term variabiliti es of the progenitors , 
as those of V838 Mon and V4332 Sgr (I Poranskii etal] 120071: 
[Kimeswenger 2007), which were sometimes raised as evidence 
against the merger hypothesis, appear now natural in view of our 
data for the V1309 Sco progenitor We do not claim that all ob- 
served eruptions of the V838 Mon-type are mergers of contact 
binaries. There can be different ways leading to stellar mergers. 
What the case of VI 309 Sco evidently shows is that the obser- 
vational appearances of a stellar merger are indeed the same as 
those observed in the V838 Mon-type eruptions. 

The outburst of VI 309 Sco was shorter and less luminous 
than those of V838 Mon and the extragalactic red novae. The 
latter objects attained luminosities of about or above 10^ Lq and 
their eruptions lasted a few months. These differences are most 
likely caused by the masses of merging stars. For V8 38 Mon, an 
~ 8 Mq primary was probably involved dTvlenda & Soker 2006) 
instead of a ~ 1 Mq one of VI 309 Sco. 

As noted in Sect. 15.2.31 the orbital period of the VI 309 Sco 
progenitor of ~1.4 day is long for the observed population of 
contact binaries. From a study of contact binaries discovered 
by the OGLE proj ect in a sky reg ion very close to the posi- 
tion of VI 309 Sco. iRucinsd (119981) concluded that the W UMa 
type seq uence sharply ends at t he orbital period of 1.3-1.5 days 
(see also iPaczvnski et ani2006h . i.e. just at the orbital period of 
the VI 309 Sco progenitor. This can be a pure coincidence (just 
one case observed), but can also indicate that binaries passing 
through contact at periods k. 1 day are not rare, but that the con- 
tact phase in their case is relatively short and quickly leads to a 
merger 

VI 309 Sco, an overlooked object (only one research paper 
published so far, i.e. Mason et al. 2010), deserves much more 
attention of the observers and astrophysicists, as do the other 
V838 Mon-type objects. Apart from supernovae, they belong to 



the most powerful stellar cataclysms. As often happens in nature, 
cataclysms destroy old worlds, but also give birth to new ones. 
What will develop from the stellar mergers? Fast rotating giants, 
similar to FK Com? Peculiar stars with circumstellar discs, when 
new generation planets can be formed? To answer these ques- 
tions, we just have to follow the evolution of the V838 Mon-type 
objects, VI 309 Sco in particular. 
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Appendix A: Periodograms from the preoutburst 
observations 

Figures lATT] and lAl2l present periodograms derived from the OGLE observations 
of the progenitor of VI 309 Sco in 2002 - 20 07. The per iodograms were ob- 
tained using the method of Schwarzenberg-Cz emvl 119961) . which fits periodic 
orthogonal polynomians to the observations and evaluates the quality of the fit 
with an analysis of variance (AOV). The resulting AOV statistics is plotted in the 
ordinate of the figures. 

In all seasons the dominating peaks is at a frequency of ~0.7 day"', which 
we interpret as a frequency of the orbital period of the contact binary progeni- 
tor (Sect. 15.2) . In earlier seasons, there is also a strong peak at a frequency of 
~1.4 day"', which is because the observations can also be reasonably fitted with 
quasi-sinusoidal variations having a period twice as short as the orbital period. In 
2006 this peak is much weaker and practically disappears in 2007, which reflects 
the evolution of the light curve as displayed in Fig. [5] All other peaks are aliases 
resulting from combinations of the two above frequencies and a frequency of 
1 day"' or subharmonics of the main peaks. 
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Fig. A.l. Periodograms from the observations of the progenitor 
of VI 309 Sco in 2002-2004. 
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Fig. A.2. Periodograms from the observations of the progenitor 
of V1309 Sco in 2005-2007. 



